Coherent regimes of mutually coupled Chua's circuits 
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We study the dynamical regimes that emerge from the strongly coupling between two Chua's 
circuits with parameters mismatch. For the region around the perfect synchronous state we show 
how to combine parameter diversity and coupling in order to robustly and precisely target a desired 
regime. This target process allows us to obtain regimes that may lie outside parameter ranges acces- 
sible for any isolated circuit. The results are obtained by following a recently developed theoretical 
technique, the order parameter expansion , and are verified both by numerical simulations and on 
electronic circuits. The theoretical results indicate that the same predictable change in the collec- 
tive dynamics can be obtained for large populations of strongly coupled circuits with parameter 
mismatches. 
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I. INTRODUCTION 

One of the most important features of interacting oscil- 
lating units is the possibility for their motions to entrain 
each other and evolve in synchronicity. The phenomenon 
of synchronization is a universal feature underlying the 
collective behavior of populations of dynamical systems 
and its importance has been recognized in many different 
fields 0,0 13- In particular, a great deal of attention has 
been devoted to globally-coupled dynamical systems as 
a model for a variety of physical Q , chemical 0| and 
biological systems 0, ■ 

Contrary to the assumption of identical units, common 
in many theoretical approaches, real systems display an 
unavoidable diversity in the oscillator population show- 
ing up as, e.g., parameter mismatches. The effect of such 
a diversity might hinder synchronization and yields an 
incoherent regime if the coupling is too weak. The main 
feature of incoherence is that the motions of the individ- 
ual systems are largely independent to each others. The 
collective dynamics, then, is either stationary in average 
(for infinite number of elements) or displays fluctuations 
that scale with the population size. When the coupling 
strength is increased, some or all the elements of the pop- 
ulation synchronize and, as a result, oscillations start to 
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be detected at a macroscopic level. The transition to 
collective oscillations is typically accompanied by rather 
complex regimes, such as clustering |p , partial synchro- 
nization [lfj, phase synchronization |l 1| . etc. For even 
stronger coupling, complete synchronous regimes, char- 
acterized by a high degree of coherency within the pop- 
ulation, develops. This means that all the population 
elements share the same kind of dynamics, apart from 
small differences due to the microscopic diversity, and 
this dynamics reflects at the macroscopic level. 

Despite of the numerous theoretical advances in the 
analysis of globally coupled dynamical systems, the ex- 
perimental verification has revealed to be extremely dif- 
ficult, in particular as far as biological systems are con- 
cerned. The Kuramoto transition from incoherence to 
the locked state is in this sense a good example It 
dates back to 1975 and is one of the best known and gen- 
eral results in synchronization theory. However, to our 
best knowledge, it has been quantitatively verified once 
in experiments and required the implementation of an 
ad hoc system |l3f . Two main problems are encountered 
in experiments. One is the difficulty in controlling or 
measuring the parameter variability in the system. The 
second one is related to the global coupling which is very 
rare in nature. These experimental limitations are espe- 
cially relevant for regimes close to the incoherence one 
where a large population is needed. 

In order to avoid these problems associated with inco- 
herence, in this paper we focus instead on globally cou- 
pled systems close to the synchronous state. As it was 
recently shown, such coherent states do not strongly de- 
pend on the population size 0, 0] . This fact is very 
interesting from the experimental viewpoint because the 
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main characteristics can be studied in systems of only two 
elements and the results extended to larger populations. 
Moreover, the coherent dynamics may be described by a 
set of low dimensional equations with only few param- 
eters. Transitions among different regimes correspond 
to bifurcations involving a small number of dimensions 
and parameters, thus having the advantage of being very 
robust to changes in the microscopic features of the pop- 
ulation. 

As an application of such theoretical results, we study 
the strong coupling regimes of two Chua's circuits. 
Chua's circuits offer an intermediate step between theory 
and applications. They can be modelled with a simple 
set of non- linear equations, thus allowing the use of ana- 
lytical methods and quick numerical simulations. At the 
same time they can be easily implemented in hardware. 
Even if they constitute a rather controlled experimental 
system, a verification of theoretical results faces a test 
on its robustness with respect to: (i) noise fluctuations 
in states and parameter values, (ii) constrains in the ac- 
cessible parameter precision and ranges, and (iii) small 
functional mismatches between model and implementa- 
tion. Finally, Chua's circuits are particularly interesting 
because they constitute prototypical devices for nonlin- 
ear system testbeds. 

Section[n]derives the equations for the macroscopic dy- 
namics of an interacting population (of any size) in the 
coherent regime. This representation allows us to infer 
on the dynamical regimes of the experimental system. 
We then provide two examples of comparison between 
theoretical prediction and experimental system. In Sec. 
IIIII we show how, by controlling parameter mismatch, 
one can drive the mean field to regimes that are different 
from the uncoupled dynamics and that might not be ac- 
cessible by any of the individual systems. In Sec. II VI we 
consider the case in which the individual elements differ 
for the time scale of their motion, and focus in particular 
on the regime of oscillaton death, where parameter diver- 
sity induces a suppression of the individual and collective 
oscillations. Finally, Sec. is devoted to the discussion 
of the results and to their possible applications. 

II. ORDER PARAMETER EXPANSION OF 
MUTUALLY COUPLED CHUA'S CIRCUITS 

The order parameter expansion [T3.[l5l | is a useful tech- 
nique for studying regimes that bifurcate from the per- 
fectly synchronous state in populations of globally and 
strongly coupled elements when the parameter mismatch 
is increased. Let us here review the main results by con- 
sidering a system of the form: 

x j = f(x j ,p j ) + K(X-x j ) j = 1,2,..., TV. (1) 

where pj is a parameter containing the diversity of the 
system and the dynamics of the each element fj(xj) = 
/ (xj ,Pj), is defined by the smooth function fj : W 1 —> 
M. n . All the elements are coupled to the the mean field of 



the population X — (xj) through the coupling function 
K(X — Xj) that, as we will see later, for the case of the 
electronic coupling we use, can be assumed as a linear 
diagonal operator (in fact, a more general system can be 
considered, see jig for details.) The basic idea of the 
method is to obtain an effective equation of motion for 
the mean field variable A, valid when all the elements 
evolve in time close to X. This is done by considering 
regimes close to the perfectly synchronous state, defined 
as Xj = X, Vj. Under this assumption, equation 
can be expanded into the individual deviations from the 
mean field tj = xj — X and from the mean parameter 
5 j = pj - p , being p Q = (p). By averaging Eq. (JTJ 
over the population, one gets, at first order, the following 
approximated equation for the mean field: 

'X = F(X,p ) + D x . p F{X,p ) W 

(2) 

W = a 2 D p F(X,p ) + {D x F(X, Po ) - K) W, 

where W = (5 e) is a second macroscopic variable, or 
shape order parameter; po is the average parameter value 
and a 2 is the standard deviation of the parameter distri- 
bution. Besides being low dimensional, such equation has 
the advantage of having a clear physical meaning: the 
mean field behaves like the average element, perturbed 
by a macroscopic variable that quantifies the mismatch 
in phase and parameter space. The macroscopic regimes 
thus appear as an unfolding of the perfectly locked state 
X = F(X,po) and are controlled by the parameter K 
and a 2 . Such equation holds when the displacements 
from the mean field ej are small, and thus, operationally, 
when the coupling is sufficiently strong compared with 
the parameter mismatch. Another important considera- 
tion is that Eq. @ does not depend on the population 
size. Consequently, the results we obtain for two Chua's 
circuits can be readily generalized to larger arrays. 

Let us now apply Eq. to the case of two mutually 
coupled Chua's circuits. The configuration of the circuits 
that we study is shown in Fig.^ The equations are easily 
written in terms of the voltages v\ and V2 (at point 1 
and 2 respectively) and the current I at point 3, the first 
circuit is described by: 

where R]j 2 are the two coupling resistors. These resis- 
tors will be used as control parameters. Other internal 
parameters, that are kept constant, are displayed in the 
caption of Fig. ^ The second circuit has an equivalent 
equation, with lower index 2 instead. Before proceeding 
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FIG. 1: Scheme of the two mutually coupled Chua's circuits. 
The parameters are C\ a = 10.2 nF ,Ct 2 = 101.8 nF, Li, a = 
18.4 mH, fli, 2 = 1562 Ohm, G?, 2 = -7.625 10~ 4 , G\ a = 
4.12 10~ 4 . 



with the analysis, it is worth mentioning that the mu- 
tual coupling Eq. J3J term can be recast in the form of 
a global coupling term, so that the correspondence with 
Eq. J2J is readily established. The dependence on the 
mean field becomes apparent by noticing that: 

{v\ v\) = (v\ + v\ - 2v\) = 2 (( Vl ) - «}) . (4) 

The coupling constants are the multiplicative factors 
l/(C x R\) and l/(C 2 R 2 k ). At this point, we recover Eq. 

by identifying Xj with the vector < wj, v 2 , Ij > and 
fj(xj) with the equation for an uncoupled Chua circuit. 
The coupling matrix is: 



K 1 



















1 



(5) 



Since we assume identical capacities for the two cir- 
cuits, the coupling is symmetrical. The explicit calcula- 
tion of the reduced equations are given in the Appendix. 
With this formalism, we are now ready to use Eq. (J2J 
on some specific condition and to study the collective 
regimes that emerge due to the interplay of parameter 
mismatch and coupling. 



III. STRONG MUTUAL COUPLING 

Let us first address the case in which the mismatch 
is in the internal resistors R\ t 2 of the two circuits, while 
the other internal parameters are kept identical. The cou- 
pling resistors R^' 2 are chosen weak enough to ensure that 



the coupling is strong and the regime coherent. Under 
such strong coupling conditions, the collective dynamics 
will be qualitatively equal to the dynamics of an "aver- 
age oscillator" . This corresponds to a system of equations 
identical to those of each single element, except for the 
parameter on which diversity is imposed, that is instead 
substituted by its average value. Intuitively, such prop- 
erty of the macroscopic dynamics is due to the fact, that 
when the coupling matrix is predominantly diagonal, the 
shape parameter W is kept small. This corresponds to 
the situation in which the population is tightly packed 
around the mean field. In the limit in which the coupling 
is infinitely strong, W vanishes and in Eq. (0) we end up 
with the equation for an average oscillator. 

Such a consideration has interesting implications in the 
control of the collective behavior of populations of dy- 
namical systems, as well as in the detection of the micro- 
scopic features in terms of purely macroscopic observa- 
tions. Indeed, the dynamics that is stabilized by coupling 
non-identical circuits may not be attainable by any of 
them if uncoupled. The collective regime is thus "coded" 
in the diversity and coupling rather than in the inter- 
nal features of the individual circuits. Moreover, Eq. @ 
allows a straightforward identification of such a hidden 
collective behavior. 

As a first application of the aforementioned theoreti- 
cal results, let us examine the case when both circuits 
operate in a periodic regime. Figure |21 (a) and (b) dis- 
plays the experimentally measured attractors of the two 
uncoupled circuits, of period 3 and 4 respectively. Figure 
|51(c) and (d) show that, as a result of the coupling, both 
oscillators behave chaotically, although the mismatch in 
the parameters has only minor effects on their individual 
dynamics. 

The same idea can be used for stabilizing unstable 
periodic orbits. Let us, for instance, consider the case 
in which both circuits are chaotic, but their parameter 
values are located, when the circuits are uncoupled, on 
the opposite side of a periodic window. Figure [3] shows 
that once they are coupled, both circuits get synchronized 
onto a period-three orbit. The coupling thus results into 
the stabilization of an unstable orbit embedded into the 
chaotic attractor. 

The possibility of obtaining a qualitatively novel dy- 
namics by coupling two circuits with different internal 
parameters is robust with respect to changes in the cou- 
pling intensity (also under differential modifications of 
the resistors) as long as both coupling resistors are suffi- 
ciently weak, that is coupling is sufficiently strong. This 
robustness was confirmed by repeating the two experi- 
ments described above for several values of the coupling 
resistors. 

Besides providing a way for targeting specific regimes, 
the experiments also show that the mathematical limit of 
infinitely strong coupling is actually attainable, to a good 
approximation, even for experimental feasible resistance 
values. It is also important to remark that in order to 
target a regime one only needs to know the bifurcation 
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FIG. 2: Collective chaotic behavior as generated by coupling 
two nonidentical circuits, each of which is individually pe- 
riodic, a) and b) Uncoupled chaotic dynamics of the two 
circuits (Rj = 1740 Ohm, R 2 = 1644 Ohm), c) Attractor of 
the first circuit after the coupling is established. Being the 
coupling strong (Rj,' 2 = 200 Ohm), the attractors of the sec- 
ond circuit and of the mean field are nearly identical to this 
one and thus are not shown, d) The same behavior of the 
strongly coupled circuits is predicted by a single circuit with 
the internal resistance set at (Ri + R-i)j1 = 1692 Ohm. 



FIG. 3: Stabilization of an unstable period-three orbit by cou- 
pling of two nonidentical chaotic circuits, a) and b) Uncou- 
pled chaotic dynamics of the two circuits (R\ = 1621, R2 = 
1745 Ohm), c) Attractor of the first circuit after the coupling 
is established. Being the coupling strong, the attractors of 
the second circuit and of the mean field are nearly identical 
to this one and thus are not shown, d) As in the FigE?lthe 
dynamics of the coupled circuits can be predicted by a sin- 
gle circuit with the internal resistance at the average value 
{R 1 + R 2 )/2 = 1683 Ohm. 



diagram: this information can be obtained experimen- 
tally and in principle does not require the knowledge of 
the functional form of /. 



IV. OSCILLATION DEATH IN CIRCUITS 
WITH TIME SCALE MISMATCH 

In this section we address oscillaton death, a regime 
where the dynamics of each population element is sup- 
pressed as a collective effect, i.e., resulting from the in- 
terplay of coupling and diversity. Oscillator death has 
been firstly described in populations of limit cycle os- 
cillators with direct 0, Hsj or time-delayed coupling 
|19| . Recently, an order parameter expansion has been 
used to demonstrate that this regime may arise under 
generic conditions in populations of globally coupled el- 
ement with time scale mismatch [l5j. Such results in- 
troduce us to the next study on coupled circuits. The 
equations for a population with time scale mismatch and 
global, linear coupling can be written as: 

Xj=t j J{x ] ) + K{X-x j ), j 1.2 V. (6) 

This equation may be seen as a simple way for gener- 
alizing the Kuramoto model, where the distributed pa- 
rameter is the natural frequency of the oscillators, to an 
individual dynamics different from a limit cycle. The pe- 
culiar parameter dependence of Eq. JJJJ is reflected in a 



simple form of the reduced system: 
(x = (r)f(X) + JW 

(7) 

{w = a 2 f(X) + ({T)J-K)W, 

where J is the Jacobian D x f(x). 

Due to the multiplicative nature of the distributed pa- 
rameter in Eq. ©, any fixed point of / is also a fixed 
point for the mean field. The stability of such a point, 
however, can change due the presence of the coupling and 
time scale mismatch. With some algebra one finds that 
generic conditions exist under which an unstable focus 
of f becomes attracting for the population. Referring to 
|15| for the details, here we just notice that, if the time 
scale mismatch is sufficiently wide, this equilibrium is 
stabilized for "intermediate" values (i.e., large enough to 
avoid incoherence and smaller than those where nonsta- 
tionary coherent solutions are present) of the coupling. 

For the case of Chua circuits, one can see from Eq. © 
that the time scale can be easily changed by tuning the 
capacitors Ci, C2 and the inductance L. We modify the 
time scale of one of the two circuits by making use of 
commercially available capacitors and inductances, thus 
obtaining a mismatch of at most 14%. Both circuits, 
when uncoupled, have similar chaotic dynamics (double 
scroll attractor) although their time evolution takes place 
with different speed. 

The bifurcation diagram of the reduced system Eq. 
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FIG. 4: Bifurcation diagram of the reduced system Eq. £J. 
The Poincare section of the attractor is shown as a function of 
the parameter mismatch between the two circuits, indicating 
that an increase in diversity induces a transition from a two- 
lobes to a one-lobe attractor. This macroscopic attractor then 
undergoes a backward period-doubling bifurcation cascade to 
a limit cycle, and eventually a fixed point is stabilized through 
a subcritical bifurcation. 

can be numerically computed for an interval of time 
scale differences comparable to that accessible in our ex- 
perimental setup. The diagram provides the coherent 
regimes of diffusively coupled Chua circuits as a function 
of their parameter diversity and is indipcndent of the 
population size. As expected from the previous section, 
for a low value of the coupling resistors (strong coupling) , 
the circuits are entrained on a chaotic attractor indistin- 
guishable from the uncoupled dynamics. However, when 
the coupling is reduced a cascade of coherent regimes bi- 
furcates from the chaotic attractor and eventually the 
state of oscillator death is reached (Fig. . 

The resistances R\—R\ — Rk of the coupled circuits 
can now be tuned in order to target a specific collective 
regime, according to the numerical characterization of 
the reduced system. Figure[S]shows some of the emergent 
dynamics experimentally observed. The oscillation death 
regime also appears in the circuits for R k slightly above 
1.5 kf2 (not shown). 



V. CONCLUSIONS 

Although the importance of emergent behavior in pop- 
ulations of nonidentical dynamical systems is widely rec- 
ognized, only few experimental systems offer the neces- 
sary degree of control for testing the theoretical results. 
Here, we have proposed the implementation of a circuit 
aimed at mimicking the coherent behavior of populations 
with parameter mismatch in the region of strong cou- 
pling. Theoretical results and numerical simulations sug- 



FIG. 5: Experimental behaviour of the coupled circuits at 
maximum coupling Rk — (a) and at values indicated in the 
bifurcation diagram of the reduced system Fig.0](b-d). 



gest that this case is representative of larger populations 
and constitute a tool for addressing diversity-dependent 
coherent dynamics. 

The experimental system was composed of two sym- 
metrically coupled Chua circuits, the diversity in their 
components introduces a mismatch in the internal pa- 
rameters. We have shown that the average dynamics of 
such coupled systems can be very different from that of 
the individual circuits when decoupled. The order pa- 
rameter expansion developed for populations of strongly 
coupled dynamical systems allowed us to predict such a 
collective behavior based on the knowledge of the uncou- 
pled dynamics. In this way, the emergent dynamics can 
be tuned, by controlling the parameter mismatch, onto 
an attractor that might not be reachable with any of the 
individual circuits. The possibility of robustly targeting 
a specific coherent regime can find a wide applicability 
in populations of globally coupled non- identical systems, 
the chaotic behavior of Chua circuits being just one ex- 
ample of individual dynamics. 

As a second example of diversity-induced qualitative 
changes of the collective behavior, we have addressed the 
case when the circuits differ by their time scales. In this 
case, the assumption that the coupling matrix is inde- 
pendent of the parameters is valid for sufficiently small 
parameter mismatch. 

The coherency of the dynamics is not only the basic 
assumption that allows us to analytically address the col- 
lective regimes, but also guarantees that the results ob- 
tained for two circuits can be extended to populations of 
many elements with a comparable variance of the param- 
eter diversity. The robustness of strong coupling regimes 
to changes in the number of interacting circuits is con- 
firmed by numerical simulations. These show indeed that 
the collective behaviour is selected according to the pa- 
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rameter distribution variance rather than from the actual The coupling matrix K, is: 
population size. 
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VII. APPENDIX 

The equations approximating the average dynamics of 
a population with time scale mismatch can be computed 
for the Chua circuits that we have addressed in this pa- 
per starting from Eq. The collective dynamics is 
described by a system of six equations (three for the vari- 
ables of the average X — [(v\), (1)2), (cr)] and three for 
the shape parameter W. This can be compactly written 
as ©, where f(X) is a vector of components: 



J 



where 



3i,i 
1 

C 2 R 



1 


" 


Ci R 


1 


1 


C 2 R 


c 2 


1 


r b 


L 


L J 



/(i)= -7^1 



/(2) = 
/(3) = 



v 2 - V\ 



G a - G b 



C 1 RC 1 2C 1 

x (\ Vl +BP\ - \v!-BP\) 

v 2 s 



- (l ' 11 

C 2 



1 
L 



Ir b 



R 



3i,i 



"ci 



Rd 



G a — Gb 
2Ci 



x [sgn(vi + B p ) - sgn(vi — B p )] 



[1] A. Pikovsky and Y. Maistrenko, eds., Synchro- 
nization: Theory and application (Kuwler, Dor- 
drecht/Boston/London, 2003). 

[2] E. Mosekilde, Y. Maistrenko, and D. Postnov, Chaotic 
Synchronization: Applications to Living Systems (World 
Scientific, Singapore, 2002). 

[3] S. C. Manrubia, A. S. Mikhailov, and D. H. Zanette, 
Emergence of dynamical order: synchronization phenom- 
ena in complex systems (World Scientific, Singapore, 
2004). 

[4] R. A. Oliva and S. H. Strogatz, Int. J. of Bif. and Chaos 
11, 2359 (2001). 

[5] Y. Kuramoto, Chemical Oscillations, Waves and Turbu- 
lence (Springer, Berlin, 1984). 

[6] I. Z. Kiss, Y. Zhai, and J. L. Hudson, Science 296, 1676 
(2002). 

[7] S. Dan0, P. G. S0rensen, and F. Hynne, Nature 402, 320 
(1999). 

[8] A. T. Winfree, The Geometry of Biological Time, 2nd 

edition (Springer, New York, 2001). 
[9] G. Osipov and M. Sushik, Phys. Rev. E 58, 7198 (1998). 



[10] P. C. Matthews, R. E. Mirollo, and S. H. Strogatz, Phys- 
ica D 52, 293 (1991). 

[11] M. G. Rosenblum, A. Pikovsky, and J. Kurths, Phys. 
Rev. Lett. 76, 1804 (1996). 

[12] Y. Kuramoto, in International Symposium on Mathemat- 
ical Problems in Theoretical Physics, edited by H. Araki 
(Springer, New York, 1975), vol. 39 of Lecture Notes in 
Physics, p. 420. 

[13] I. Z. Kiss, Y. Zhai, and J. L. Hudson, Science 296, 1676 
(2002). 

[14] S. De Monte and F. d'Ovidio, Europhys. Lett. 58, 21 
(2002). 

[15] S. De Monte, F. d'Ovidio, and E. Mosekilde, Phys. Rev. 

Lett. 90, 054102 (2003). 
[16] S. De Monte, F. d'Ovidio, H. Chate, and E. Mosekide, 

Phys. Rev. Lett. 92, 254101 (2004). 
[17] G. B. Ermentrout and W. C. Troy, SIAM J. Math. Anal. 

20, 1436 (1989). 
[18] G. B. Ermentrout, Physica D 41, 219 (1990). 
[19] F. Atay, Phys. Rev. Lett. 91, 094101 (2003). 



